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I hereby declare that the nonprofit organization identified above qualifies as a nonprofit organization as defined in 37 
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to the invention entitled METHODS FOR ESTIMATING THE BODY VOLTAGE OF D IGITAL PARTtAf.f.Y 
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NAME 

ADDRESS 

[ ] INDIVIDUAL [ ] SMALL BUSINESS CONCERN [ ] NONPROFIT ORGANIZATION 

NAME 

ADDRESS 

[ ] INDIVIDUAL [ ] SMALL BUSINESS CONCERN [ ] NONPROFIT ORGANIZATION 

I acknowledge the duty to file, in this application or patent, notification of any change in status resulting in loss of 
entj|ement to small entity status prior to paying, or at the time of paying, the earliest of the issue fee or any 
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BAKER BOTTS, LX.P. 
30 ROCKEFELLER PLAZA 
NEW YORK, NEW YORK 101 12 



TO ALL WHOM IT MAY CONCERN: 

Be it known that I, Kenneth Shepard, a citizen of the United States of America, 
residing in City of Ossining, County of Weschester, whose post office address is 46 Snowden 
Avenue, Ossining, NY 10562, respectively, have invented 

METHODS FOR ESTIMATING THE BODY VOLTAGE OF DIGITAL 
PARTIALLY DEPLETED SILICON-ON-INSULATOR CIRCUITS 

of which the following is a 

SPECIFICATION 

NOTICE OF GOVERNMENT RIGHTS 
The U.S. Government has certain rights in this invention pursuant to the terms of the 
National Science Foundation award number CCR-97-34216. 

BACKGROUND OF THE INVENTION 
I. Field of the invention. 

The present invention relates to partially depleted silicon-on-insulator technology, 
and more particularly to techniques for estimating the body voltage of digital integrated circuits and 
partially depleted silicon-on-insulator ("PD-SOI") technology. 
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IL Description of the related art. 

Silicon-on-insulator ("SOI") technology has long been used for radiation-hardened 
or high- voltage integrated circuits. In more recent years, SOI has also been considered for us as a 
technology for high-performance, low-power deep-submicron digital integrated circuits. For digital 
applications, fully-depleted devices have been largely abandoned in favor of partially-depleted 
technology, because of the difficulty in controlling the threshold voltage of fully-depleted thin-film 
transistors. Partially-depleted SOI ("PD-SOI") technology has two main advantages for digital 
applications: the reduction of parasitic source-drain depletion capacitances and the reduction of the 
reverse-body effect in stack structures and pass-transistor logic. 

At the device and circuit level, however, the floating body effect in PD-SOI has posed 
and continues to pose major challenges in the successful use of this technology. As discussed in 
Pong-Fei Lu et al., "Floating-body Effects in Partially Depleted SOI CMOS circuits, 32 IEEE J. 
Solid-State Circuits 1241 (1997), there is a parasitic bipolar effect which can result in noise failures 
if not correctly considered. In addition, there can be large "uncertainties" in the body potential, and 
consequently the threshold voltage, of devices due to unknown past switching activity. For many 
circuits, the design margining required to protect against this uncertainty erodes all of the potential 
performance advantage under nominal operation. In addition, for many circuit styles in which noise 
margin is strongly determined by threshold voltage (e.g. dynamic circuits), considerable overdesign 
for noise can also result from conservative body-voltage margining. 

Several circuit design techniques, such as those described in D.H. Allen et al., "A 
0.20 \im 1 .8 V SOI 550 MhZ 64b Power PC Microprocessor with Cu Interconnects," Digest Tech. 
Papers, ISSCC, pp. 438-39 (1999), and C.T. Chuang et al., "SOI Digital CMOS VLSI - A Design 
Perspective," 36th ACM/IEEE Design Automation Concerence, pp. 709-714 (1999), attempt to 
contain the noise impact of the parasitic bipolar current and the delay and noise-margin variation due 
to the floating body. Some of these design techniques, such as predischarging internal nodes of the 
nFET pulldown stack in domino logic to avoid parasitic bipolar currents, are quite counter to design 
practice in bulk. 
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In addition, previous circuit-level modeling work on PD-SOI has focused either on 
device issues, or on delay and noise effects due to the floating-body effect evident for particular 
circuits under periodic stimulus, such as pulse stretching, frequency-dependent delay time. See, for 
example, J. Gautier et al,. "On The Transient Operation of Partially Depleted SOI NMOSFET;s," 
5 161EE Electron Device Letters 498 (1 995), and R. Puri et al, "Hysteresis Effect in Pass-Transistor- 
Based Partially-Depleted SOI CMOS Cirsuits," Proc. Int'l SOI Conf. 1 998. However, neither such 
modeling work nor the circuit design techniques discussed above are able to accurately estimate the 
upper and lower bounds of the body voltage of PD-SOI circuits in order to characterize the important 
electrical characteristics of the circuit. Accordingly, there exists a need in the field to accurately 
1 0 estimate the upper and lower bounds of the body voltage of PD-SOI circuits, taking into account the 
5 past history of the circuit, so that certain important electrical characteristics, such as delay and noise 
W margins, of the circuit can be properly considered without requiring excessively long simulatrons. 

W SUMMARY OF THE INVENTION 

I §3 An object of the present invention is to provide a technique for accurately accurately 

O estimate the upper and lower bounds of the body voltage of PD-SOI circuits. 
[Lj A further obj ect of the present invention is to provide a technique for determining the 

2 body voltage of PD-SOI circuits which accounts for the past history of the circuit. 
P Yet a further object of the present invention is to provide techniques for accurately 

20 determining important electrical characteristics,such as delay and noise margins, of PD-SOI circuits. 

In order to achieve these objectives as well as others that will become apparent with 
reference to the following specification, the present invention provides a method for statically 
estimating a body voltage of one or more transistors which form digital partially depleted silicon-on- 
insulator circuit. In one technique, one or more device models are obtained and abstracted to 
25 generate one or more simplified electrical descriptions, each corresponding to one of the transistors. 
The device models may be n or p channel Field Effect Transistor models. As part of this abstraction, 
the circuit topology is then checked to generate one or more sets of accessible states, with each set 
corresponding to one of the transistors and being indicative of whether under any allowable 
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switching activity, the source, gate, or drain could be high or low. Next, sets of reference state body 
voltage minima and reference state body voltage maxima are determined for each of the transistors 
based on corresponding simplified electrical descriptions and corresponding sets of accessible states . 
Finally, a target state body voltage minima and a target state body voltage maxima are ascertained, 
5 one for each transistor, based on the determined sets of reference state body voltage minima and 
reference state body voltage maxima. This target state reflects the appropriate initial condition for 
the electrical simulation in question. The electrical simulation in question also determines whether 
the minima or maxima is used for a given FET. 

Preferably, the abstracting step includes abstracting each device model to obtain 
10 displacement voltages d i ? steady-state reference voltages Vf ero , and forward biased reference 
yj voltages V/ orward , associated with each transistor state i for each such model. Using such abstracting, 
m the determining step advantageously involves determining the minima and maxima reference state 
^ body voltage (V B ref ) m i n ° r (V B ref )m ax f° r eac h transistor in the given circuit. With such a methodology, 
HJ the present invention advantageously provides for ascertaining one or more target state body voltage 
IS minima, and one or more target state body voltage maxima, based on the determined reference state 
Jf? body voltage minima and maxima using either foil uncertainty or accessibility analysis, and the 
□ displacement abstractions d r 

q In an especially preferred technique, an additional step of performing active net 

tagging on each net in the circuit to determine whether any of nets will switch with regular frequency 
20 is performed prior to the determination of reference state body voltage minima and maxima. In this 

arrangement, the simplified electrical abstractions are thedz K^and V/ orward combined with the 

results of active net tagging, and the determining step advantageously involves determining 
minimum and maximum reference state body voltages for each transistor in the circuit based on 
corresponding displacement voltages d i5 corresponding sets of accessible states. With this 
25 methodology, the present invention also calls for ascertaining one or more minimum target state 
body voltage minima and one or more target state body voltage maxima, based on the determined 
reference state body voltage minima and maxima using modified accessibility analysis. 
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In another preferred technique, the device models are also abstracted to obtain time 
constant characterizations, and an additional step of calculating signal probabilities and timing 
windows by Boolean and timing analysis is performed prior to the determination of reference state 
body voltage minima and maxima. In this arrangement, the determining step advantageously 
5 involves determining the sets of reference state body voltage minima and maxima based on 
corresponding displacement voltages d,, the time constraints, corresponding sets of accessible states, 
and the results of the calculated signal probabilities and timing windows. With such a methodology, 
the present invention calls for ascertaining the target state body voltage minima and maxima, based 
on the determined reference state body voltage minima and maxima using probabilistic analysis. 
10 In order to take advantage of the ability to accurately determine target state body 

O voltage minima and maxima, the present invention provides techniques for analyzing an electrical 
LR property of a digital partially depleted silicon-on-insulator circuit. In one preferred methodology, 
jjf the technique involves ascertaining a target state body voltage minimum and a target state body 
W voltage maximum for each transistor in the circuit, establishing an initial condition for the circuit 
1S3 by selecting either the target state body voltage minimum or the target state body voltage maximum 
O for each transistor, applying a voltage to the circuit, and measuring the sought after electrical 
'tl property of the circuit. 

'h In accordance with the present invention, the target state body voltage minima and 

S maxima are ascertained alternatively from full uncertainty estimation, accessibility analysis, 

20 modified accessibility analysis, or probabilistic analysis. If the sought after electrical property is a 
switching delay, the measuring step involves measuring a delay between a switching input and a 
switching output. Such delay measurements can be used as constituent simulations in static timing 
analysis. If the sought after electrical property is noise, the measuring step involves measuring noise 
on one or more nets in the circuit such noise measurements can be used as constituent simulations 

25 in static noise analysis. In its ability to accurately determine target state body voltage minima and 
maxima, the present invention fulfills the need to accurately estimate the upper and lower bounds 
of the body voltage of PD-SOI circuits, taking into account the past history of the circuit, so that 
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certain important electrical characteristics, such as delay, induced noise, and noise margins, of the 
circuit can be properly considered. 

The accompanying drawings, which are incorporated and constitute part of this 
disclosure, illustrate a preferred embodiment of the invention and serve to explain the principles of 
the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a state diagram for a partially depleted silicon-on-transistor pFET; 
Figs. 2a and 2b are graphs displacements for an example technology; 
Fig. 3 is an illustrative drawing showing the application of modified accessibility 
analysis to an exemplary circuit; 

Fig. 4 is an illustrative diagram showing an average cycle divided into a series of time 

slices; 

Figs. 5a-c are graphs plotting the relaxation times r } for an nFET with a 2.5-V supply 
as a function of local body voltage; 

Figs. 6a and 6b are graphs plotting reference body voltage against effective time; 
Fig. 7 is a flow diagram of the preferred process implemented by the present 

invention; 

Fig. 8 is an illustrative diagram showing a chain of inverters; 
Figs. 9a and 9b are illustrative diagram showing an exemplary circuit; and 
Figs. 10a - 13c are graphs illustrating the results of a comparison between the 
techniques described herein and a simulation for the an exemplary circuit. 
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DESCRIPTION OF PREFERRED EMBODIMENTS 
In order to better understand the embodiments disclosed herein, the PD-SOI device 
physics which underlie the present invention will first be explained. It should be noted that the 
following description is by way of example and not by limitation, with certain embodiments being 
described in order to best explain the principles of the present invention. 

The body potential of a PD-SOI Field-Effect-Transistor ("FET") is determined by 
capacitive coupling of the body to the gate, source, and drain, by diode currents at the source-body 
and drain-body junctions and by impact ionization, currents produced by current flow through the 
device (referred to as the on-state impact ionization current). The impact ionization currents have 
a strong supply-voltage dependence, decreasing with decreasing supply voltages and going to zero 

below 1.1 -V supplies. 

It is convenient to distinguish "fast" and "slow" processes. Fast processes can 
change the body potential on time scales on the order of or less than the cycle time, while slow 
processes require time scales much longer than the cycle time (up to milliseconds) to affect the body 
voltage. There are two fast mechanisms at work: switching transitions on the gate, source, or drain 
which are capacitively coupled to the body (which will hereinafter be referred to as coupling 
displacements), and forward-bias diode currents across source-body and drain-body junctions with 
voltages exceeding the diode turn-on voltage (which will hereinafter be referred to as body 
discharge). The slow processes involve charging or discharging the body through reverse-biased 
or very weakly forward-biased diode junctions and through impact ionization. 

Figure 1 is a state diagram abstraction for an nFET which is useful to model the 
switching history determining the body voltage of a particular device. As a usually dynamic circuit 
node, the floating body has "memory." The states 110, 120, 130, 140, 151, 1 52 denoted with solid 
circles represent "static" states, states in which the FET can be stable. The states 161, 162 are 
"dynamic" states, which are only present transiently during switching events. For example, state 1 1 0 
corresponds to the case in which the gate is high and both the source and drain are low. Arrows 
indicate possible state transitions produced by switching events in the circuits containing these FETs. 
These switching events can represent transitions from the logic state at the end of the previous cycle 
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to the logic state at the end of the current cycle or can represent hazards that occur transiently within 
a cycle. States 151 and 152 can usually be treated equivalently as the same state 150; similarly 
states 161 and 162 can usually be treated equivalently as the same state 160. 

While Figure 1 applies to an nFET, is should be noted that the state diagram of a 
pFET is similar, in which the gate is high rather than low in states 1 30, 1 40, and 1 50; and low rather 
than high in states 110, 120, and 160. 

If the device is allowed to remain in one state for a very long time, the body voltage 
in each state will achieve a dc value, hereinafter referred to as si. The dc voltages in states 1 10 and 
130 (s, and s 3 ) are zero, while the dc voltages in stages 120 and 140 (s 2 and s 4 ) are given by the 
supply voltage. s 5 is determined by the steady-state balance between a weakly forward-biased 
junction drawing current from the body and a reverse-biased junction, leaking current to the body, 
enhanced by GIDL currents. Similarly, s 6 is determined by the steady-state balance between a 
weakly forward-biased junction drawing current from the body and charging current due to reverse 
leakage of the other diode junction and on-state impact ionization. 

As an example, commercially available BSIM3SOI models for an IBM partially- 
depleted SOI technology, having a 0.25 pan effective channel length, 5 -ran gate oxide, 350 - ran 
back oxide, and 140 - nm thin silicon film are considered herein. Those skilled in the art will 
understand that the present invention applies with equal force to other commercially available and 
custom fabricated partially-depleted SOI technology. 

Using exemplary supply voltages of - 1 .0V and 2.5V, Table I lists the values of s, 
for each of the six states. The remainder of the information listed in Table I will be explained below. 
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TABLE I 






nfet 




pfet 




nfet 






pfet 








si 




forward 

V i 


si 


v zero 




si 


T zero 
v i 


forward 


si 


rrZero 
y i 


forward 

V 


l 


0 


1.77 


2.41 


0.0 


0.0 


-0.6 


0.0 


0.53 


1.13 


0 


0 


-0.6 


2 


2.5 


2.50 


3.1 


2.5 


0.73 


0.12 


1.0 


1.0 


1.6 


1.0 


0.469 


-0.127 


3 


0 


2.15 


2.85 


0.0 


-0.92 


-0.98 


0 


0.87 


1.5 


0 


-0.75 


-0.907 


4 


2.5 


3.43 


3.5 


2.5 


0.35 


-0.30 


1.0 


1.75 


1.906 


1.0 


0.14 


-0.5 


5 


0.437 


2.58 


2.79 


2.068 


-0.045 


-0.23 


0.278 


1.12 


1.47 


0.722 


-0.118 


-0.47 


6 


0.711 


1.60 


1.60 


1.976 


1.24 


1.24 


0.357 


0.61 


0.85 


0.722 


0.467 


0.17 



C In the absence of body discharge, the coupling displacements that occur with each 

5 transition in Figure 1 are completely reversible on "fast" time scales; that is, if one begins in state 

(J 1 1 0 and traverses the state diagram, returning to state 1 10 on a time scale much faster than any of 

^ the "slow" leakage mechanisms, the body voltage on return will be the same as the initial body 

O voltage, a simple result of charge conservation. Because of this, the charge stored on the body can 

O be represented as the value of the body voltage in one particular state of Figure 1 (hereinafter 

J referred to as the "reference" state), which is chosen to be state 1 20 for the nFET and state 1 1 0 for 

2| the pFET. 

O From this reference body voltage {y r B ef ), we can then determine the corresponding 
body voltage in each state i \V l B ) according to the relationship: 



V B =V B f+d i( V B f ) (1) 

25 From (1), The displacements, d ^ V s f ) ? are explicitly shown to be dependent onthe reference body 
voltage because of the strong voltage-dependence of the capacitances of the source-body, drain- 
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body, and gate-body, d, is independent of device width W because the dominant components of 
capacitance scale proportionately with W. 

Referring next to Figures 2a and b, the displacements are shown as a function of 

V B f for the example technology at the 2.5-V supply, with Figure 2a illustrating an nFET and Figure 

5 2b illustrating a pFET. In many cases, these values are not easily determined with direct simulation 
of a transition from the reference state. For example, reference body voltages for the nFET greater 
than about 3.1V would result in strongly forward-biased body- drain and body-source junctions 
making it difficult to distinguish the displacement due to the transition from state 1 20 from the body 
discharge. As a result, these curves are instead determined by isolated evaluation of the MOS C-V 

10=s model. 

IKS.* 

£ With the reference body voltage as a "state-independent" way of representing the 

* charge trapped on the body, each state / in Figure 1 may be characterized by two values of this 

[V reference voltage, V. zer and V j forward . Referring again to Table I, V™° represents the steady-state 

3 value of the reference body voltage, achieved by remaining in state i for a long time. This follows 

1 S3 immediately from the s, in each state: 

~~ -rrzero . i ft r zero \ (2) 

□ V i forward represents the value of the reference body voltage for the nFET (pFET) to which the body 

would be very quickly pulled down (up) as a result of body discharge (charge), if state / were 

accessed with a higher (lower) reference body voltage than v/ orward . The values shown in Table I 

20 presume that the fast body discharge will bring the forward-biased- junction bias down to a turn-on 
voltage of 0.6 V; it is also important to note that fast body discharge can trigger parasitic bipolar 
leakage between source and drain for FETs in state 150. This means, for example, that if a FET 

which reached a dc steady-state in state 140 (with a V r B ef of 3.3 V) switches into state 120, the 
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reference body voltage will quickly discharge to V 2 forward = 3.1 V. If the FET subsequently 
remains in state 120 for a long time, V r B ef will eventually decrease to F 2 zew = 2.5V. 

An important qualitative difference in the values between V* ero and V. forward values 

the 2.5-V and 1.0-V cases is the 6 state. At a supply voltage of 2.5 V, yf° rward = V" ro (which 

5 is not the case at the 1 .0-V supply) because of the dominating effect of the on-state impact ionization 
current. This current is so large that a strongly forward-biased junction is required to balance it in 
steady-state. 

In accordance with the present invention, once the d ; curves and the 

V zero and V forward values for each state are obtained, the target state body voltage minima and 

1Q0 target state body voltage maxima are estimated. Four illustrative estimation techniques, full 
J uncertainty estimation, accessibility analysis, modified accessibility analysis, and probabilistic 

y analysis, may be used to ascertain approximate values of the target state body voltage minima and 

O 

s target state body voltage maxima. Each is discussed in turn. 

55 In the first technique, full-uncertainty estimation, no knowledge of the switching 

1 § activity of the circuit is assumed. The maximum and minimum possible values of the body voltage 
O that cover all possible stimulus and history are chosen. A state is accessible if the circuit topology 
*** allows the state to be visited. As an example, for the nFET of an inverter, those states with the 
source high would not be accessible, because the source of the nFET is tied to ground. If A 
represents the set of all such accessible states, including possibly the dynamic state 160, the 
20 minimum and maximum body voltages are given by (3) and (4): 

r yzero (3) 

<v re f \ = max J 
V B J max j e A 



. yzero 
(y ref = min K j (4) 
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For the second technique, accessibility analysis, it is assumed that the circuit is under 
steady switching activity such that every accessible state is visited with reasonable frequency, i.e., 
on a time scale that is faster than the "slow" body-voltage mechanism. Under such conditions, the 
y forward va i ues f or t he nFET (pFET) will cap the maximum (minimum) possible value of the body 

voltage, and values for the pFET will cap the minimum possible value of the body voltage. For the 
nFET, 

mref^ = minv zero (5) 
^ V B ; min jeA j 

/vrefi =minf max v zero.mm yforward (6) 



while for the pFET, 



{V 



rej\ - maw Mn v zero:mwi v forward (7) 



static 



(V r f) =max \ : zer0 



b /j (8) 

1 s " The assumption here is that it does not matter how long an accessible state is visited; 

U it will be long enough to discharge the body down to the forward-bias turn-on voltage of the source- 
O body or drain-body diodes. In general, state 1 60 is omitted from this "accessibility" analysis ( A static 
O is used in the above equations) because state 160 is a switching state and ? therefore, cannot be 
^ assured to meet this criterion. In those cases in which the relaxation in state 1 60 is extraordinarily 
20 "fast" because of a large on-state impact ionization current, state 160 can be safely included in the 
accessibility analysis. 

While accessibility analysis does not require detailed switching knowledge, it does 
require that there is enough switching activity that every accessible state is visited with a minimum 
frequency. At times, this, too, may be difficult to ensure. The third technique, which shall be 
25 referred to as modified accessibility analysis, permits the marking of signals in the design as "active" 
signals, i.e., signals that are assured to switch with regular frequency. As those skilled in the art will 
appreciate, clock nets can usually be tagged. This active net tagging establishes a set of constraints 
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that must be satisfied by a modified accessibility set of states, A ma . These are states that must be 
visited as a result of the "active" net switches. The body voltage is determined by equations 
identical to those used for ordinary accessibility analysis, except that A statlc is replaced by A ma . Thus, 
for the nFET, 

5 

(V£ ef ) . - min V? ero (9) 
js A 

(V ref ) = min (maxV zero );min wfo^rd (10) 
B max j J 

j e A J e A 

ma 

10 while for pFET. 



1© 



(V ref ) . =max(min V z . er °; max vf ormula ) 
v B y min . . J . . J 

I 6 A / e A 



(11) 



ma 



(V ref ) = max V* ero n n\ 
v B max , j K iZ -J 



A ma will be the accessibility set for the nFET (pFET) which satisfies the constraints while ensuring 
J the maximum (minimum) value of the body voltage. As in accessibility analysis, modified 
accessibility analysis explicitly excludes the dynamic state 160 from A ma . These constraints are 

o 

□ derived from the following four rules: 

20 (1) If the gate of the FET is an active net, then there must be a state in A ma with 

the gate high and a state in A ma with the gate low. 

(2) If the source of the FET is connected to ground through a succession of 

transistors all of whose gates are active nets, then A ma must contain a state with the source low. If 

connected to supply, then A ma must contain a state with the source high. 
25 (3) If the drain of the FET is connected to ground through a succession of 

transistors all of whose gates are active nets, then A ma must contain a state with the drain low. If 

connected to supply, then A ma must contain a state with the source high. 
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(4) Same-signal correlations must be considered and can result in stronger 
constraints. For example, if gate of the current FET is marked as an active net, but must be high for 
a path to ground from the source to be present through a path of active-gate FETs, then A ma must 
contain a state in which both the gate is high and the source is low. 
5 Referring to Fig. 3, modified accessibility analysis is illustrated with reference to an 

exemplary circuit. Nets C and B (marked with the arrows) are active nets. For transistor Ml, 
therefore, there must be a state in with the source and gate both high, with the gate low, and with 
the drain high. 



10 


TABLE II 




i 


si 


nfet 

rr zero 
i 


y forward 


si 


pfet 

yzero 


y forward 




1 


0 


1.78 


2.37 


0.0 


0.0 


-0.6 


if; 


2 


2.5 


2.50 


3.1 


2.5 


0.736 


0.12 




3 


0 


2.16 


2.79 


0 


-0.9315 


-0.9505 




4 


2.5 


3.43 


3.45 


2.5 


0.349 


-0.225 




5 


0.421 


2.54 


2.72 


2.079 


-0.026 


0.225 




6 


0.560 


1.3 


1.36 


2.078 


1.373 


1.15 



2 fc 

** Table II lists values of V; Zero and V ; forward for the nFET and pFET in the circuit of Fig. 



3. To find A ma that satisfies these constraints and achieves the maximum value of V r B e , refer to 

Table II and begin with the state of largest v/ orward . Adding state 140 to A ma satisfies the drain-high 
constraint and the gate-low constraint. State 120 satisfies the constraint that there must be a state 
25 with the source and drain both high. Therefore, A ma = {140, 120} and (^ e/ ) max = 3.1 V. 

The fourth and most complex technique presented herein for estimating the target 
state body voltage minima and target state body voltage maxima shall be referred to as probabilistic 
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analysis. In this technique, the estimates provided by the simple uncertainty modes are tightened 
with stochastic techniques and timing information from static timing analysis. Thus, the behavior 
of the body over a long period of time is determined by reference to an average cycle repeated over 
and over. 

In probabilistic analysis, two average cycles are characterized, one to minimize and 
one to maximize the body voltage. Referring to Figure 4, such an average cycle is shown, divided 

into a series of time slices tf , which characterize the amount of time per cycle on the average that 

the FET spends in state /. The time slices sum to the cycle time, in accordance with equation 13: 

6 < 13) 

7 t eff -t 

i=l i cycle 

In each time slice /, the reference body voltage is relaxing to the target value 
V zer0 with a characteristic relaxation time denoted by T r The body-voltage-dependent time 

constants for this relaxation are denoted by the t } and can be captured from the device models as part 
of a technology precharacterizationm, which may be a piecewise-linear representation of the 

logarithm of t } as a function of y™ f , since varies over several orders of magnitude. Moreover, 

zj is independent of device width W because both the body currents and body capacitance scale 
proportionately with W. 

Except for the fast discharge associated with a source-body or drain-body junction 
that becomes strongly forward biased, these time constants are much larger than t cycle , and any 
voltage change during a single time slice would be imperceptible in Figure 4. From the foregoing, 
the reference body voltage at the end of the cycle V n+1 can be related to the reference body voltage 
at the beginning of the cycle V n by equation 14: 

, eff/t. 6 7prn -t.eff/t: 6 eff/t. ( 14 ) 
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The steady-state solution of this difference equation (in the approximation that the r, are much 
greater than t^) is given by equation 15: 

v ref ' - 1 » I (15) 

° If , tff/Ti 



Because of the V r B ef -dependence of the x„ these equations must be solved self-consistently. 

Figures 5a-c plot the relaxation times r, for the nFET at the 2.5-V supply as a 
function of the "local" body voltage (that is, the body voltage in state i rather than the reference body 
1 q_ voltage). The relaxation times x , and t 3 , for states 1 1 0 and 1 3 0, peak around a zero body voltage 
2 and then decrease for negative body voltage as the source-body and drain-body leakage currents 
;i . increase. The decrease for positive body voltage is even more substantial as the junctions become 
forward biased and draw significant current. The same trends are observed for the relaxation times 
SI t 2 and T 4> for states 120 and 140 > where the relaxation times P eak around the supply voltage, 
1 f* indicating that the source and drain are both at the supply voltage for states 120 and 140. 
S The relaxation time x 5 for state 1 50 peaks around zero. For voltages above the peak, 

Q t 5 decreases strongly as the body-source junction becomes forward biased. For voltages below the 
O peak, x 5 also decreases as both the body-source and body-drain junctions become more reverse- 
U biased. x 5 peaks at a much smaller time-constant value than x l through x 4 because even at a body 
20 voltage near the peak, one of the junctions (the drain-body junction) is strongly reversed biased and, 
therefore, leaking considerable current. 

For state 160, the time constants are considerably smaller (factor of 1000) than they 
are in state 150. This is due to the contribution of the on-state impact ionization current to the 
leakage currents from the floating body; a very low relaxation time will mean that devices that have 
25 any time in state 1 60 in their switching transients will have their body voltages effectively "pinned" 
at s 6 . At 1 .0 V, by contrast, t 5 and t 6 are comparable in magnitude. 
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In probabilistic analysis, the maximum and minimum possible V^f values from 
Equation 15 are ascertained upon consideration of all of the allowable t* . If the Tj in Equation 
15 were constant and equal, then to maximize V r B ef one would simply maximize the time in those 
states with higher V. ero (subject to constraints on the maximum and minimum possible values of 
5 the tf imposed by the stochastic analysis); similarly, to minimize V r B ef one must maximize the time 
in those states with lower V* ero . The x, however, are not constant, but because they have the 
property that they decrease rapidly for V r / different from V™ ro states with the largest 
continue to have the biggest "pull" and the V- er0 "ordering rule" continues to hold. 

U To formalize this, the rank of the state i is defined as an integer r, indicating the priority of states to 
l(lf{ maximize the body voltage. Higher ri means higher priority. For example, from Table I, for the 

P nFET at 2.5 V, r 4 = 6, r 5 = 5, r 2 = 4, r 3 = 3, r } = 2, and r 6 = 1. A complementary rank, r, , can also be 

S defined for minimizing V^ ef ; these priorities are simply the reverse of the maximum case. Such 

y ordering, which favors states with high source or drain voltages to maximize the body voltages while 
5 favoring states with low source or drain voltages to minimize the body voltages, is technology- 
1 r independent except for the relative positions of state 1 50 and 1 60, which depend on both technology 
and supply voltage. 

Referring next to Figures 6a-b, cycle times for the nFET at the 1 . 0-V supply is shown 
as an example of how the steady-state body voltage depends on the tf . The cycle time t cycle is 1 0 
nsec. Figure 6(a), assumes a hypothetical case in which only states 1 1 0 and 140 are accessible, and 
20 plots the reference body voltage as a function of tf given that tf - Ucu - tf . 



yref _ y. 
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The values of V™ ro ; v/ orward ; V 4 zero are also noted. Given the accessibility of states 

110 and 140, the range from V* er0 to v/ orward (denoted with the arrow), would be given by 
accessibility body-voltage estimation. The probabilistic body voltage estimate becomes 
Fj zer0 as tf -> 0 and monotonicallyincreases with increasing tf (r* > n). As tf -> ^(and 

fff 0 ). rapidly increases to F 4 zero . V r / can only increase above V l forward , with =0 

on the scale of the "fast" forward-biased discharge of the body. 

In Figure 6(b), the case in which only states 1 10 and 120 are accessible, is illustrated. 

V r B ef increases monotonically as tf increases (r 2 > rj. In this case, V™ 0 < V/ orward so that 
Vi ef astf -> 0 does not actually reach the upper bound of accessibility analysis. 

B 2 

To determine allowable values of tf for probabilistic analysis, information about 

the circuit environment of each transistor in the circuit, both logical and temporal, is required. The 
logical environment of each FET is therefore characterized by a set of signal probabilities which 
determine the possible states of the source, gate, and drain of the transistor at the end of a cycle. For 
the nFET, it is preferred that seven signal probabilities are used, as follows: 

(1 ) P(G) = probability that at the end of a cycle the gate is driven high; 

(2) P(DG)= probability that at the end of a cycle the drain is driven high 

given that the gate is low; 

(3) P(D G)= probability that at the end of a cycle the drain is driven low given 

that the gate is low; 

(4) P(SG) = probability that at the end of a cycle the source is driven high 
given that the gate is low; 
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(5) P(S G) - probability that at the end of a cycle the source is driven low 

given that the gate is low; 

(6) P(S|G) = P(D|G) = probability that at the end of a cycle the drain (and source) 

are driven high given that the gate is high; and 

(7) P(S\G) = P( D\G) = probability that at the end of a cycle that the drain 
(and source) are driven low given that the gate is high. 

The corresponding probabilities for the pFET are: P(G), P(S|G), P(S\G), P(D|G), 

P(b\G), P(S|G) = P(DG), P{SG)= P(DG). The term "driven low" means that there 

is a path to ground; "driven high" means that there is a path to supply. The following conditions on 
the conditional signal probabilities must hold: 

P(D|G) + P(D|G) < 1 

P(S|G) + P(S|G) <1 

P(D|G) + P(D|G) <1 O 6 ) 

P(S|G) + P(S|G) <1 

None of these have to sum to precisely 1 because of the possibility of floating nodes. 
With these signal probabilities, one can view the state diagram in Figure 1 as representing a Markov 
process with six-by-six transition matrix A. Assigning different transition matrices to the minimum 
and maximum cases yields difference equations 17 and 18, where the six-by-six matrices 
A max and A min are given 19 and 20: 



(F i } *«- max f z Jk (17) 

(18) 



( pmm ) K+l = A . (P min A 
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(19) 



, /rt max „ max „ max ~ max n max n max , ~, 

max = rP 1 P 2 P 3 P 4 P 5a P 5b ) T 
max jcfllllllj 



A . = r pmm p min p mm p mm p mm p mm ;r 

min xrmiii; (20) 



10 In equations 17-20, (P)k is the probability of being in state * at the end of cycle fc, 

and p } is the probability of making a transition to state i and follows directly from the source, gate, 



max 



and drain signal probabilities. For example, for the nFET, p 3 =(1- P(G))P(S G)P(DG) 



mm 



Jf while P = (1- P(G))(1- P(SG))(l- P(DG)). The maximum case assumes that the 

O floating node condition on the source or drain takes a high voltage value; the minimum case assumes 
1 fg that the floating node condition on the source or drain takes a low voltage value. 
S Diagonalizing the matrix A and finding the eigenvector associated with eigenvalue 

^ 1 (normalized so that the sum of the elements of the vector is 1) gives the steady-state values of the 
From these probabilities, one can calculate a set of thirty-six transition probabilities (that is, the 
probability that at given FET in a given cycle is transitioning from state i to state j) for both the 
20 minimum or maximum cases, in accordance with equation 21 : 

max max max min min min 

P. =p . p and P = i> P . (21) 
* -> J ' J * -> J i J J 

Next, the fraction of the cycle time { j^ k , that can be spent in each state i is 
determined as part of the transition j^k to maximize (or minimize) the body voltage among the set 
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of possible waveforms. To do this, early and late arrival times (rising and falling) for the source, 
gate, and drain of the FET under consideration (the target FET) are required. The seven early times 
are and are denoted as follows: 

( 1 ) S early . Earliest time the source of the FET can be driven high; 

V / rise 

(2) Sff . Earliest time the source of the FET can be driven low; 

(3) D e " r s l J . Earliest time the drain of the FET can be driven high; 

(4) D e ffl . Earliest time the drain of the FET can be driven low; 

(5) G°jw ■ Earliest time the gate of the FET can be driven low; 

(6) G e r ? s r J y • Earliest time the gate of the FET can be driven high; and 

(7) G"° e rly . Earliest time the gate of the FET can be driven high. 

There are comparable arrival times associated with the late state: 

S l7e > S faU ' D lTe > D fall ' G fa'll > aVld G rTe • F ° r ttle S0UrCe ^ dlain ^ FET 

is assumed to be off. 

Each transition has associated with it a set of arrival times necessary to make that 
transition, a transition set, denoted as T,_ s . For example, for the 1 10-120 transition for the nFET, 
the associated arrival time set is T^ 2 = {S rise> D n J. For 1 10- 140 for the nFET, the associated arrival 
time set is T,^ ={G fall , S nse , D ri J. The max and min operators which act on the transition sets are 
defined such that max(T H ) returns the largest of the early arrival times in the transition set, while 
min(T H ) returns the smallest of the late arrival times in the transition set. 

To indicate the states of Figure 1 involved in a cycle and to handle the possibility of 
hazards, the waveform in a cycle (in this case involving a transition from i toy) is denoted using a 
transition notation as follows: 

max 7 min 

i > k > j 
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In this cycle, a hazard to state k occurs as part of the transition. The transition notation must involve 
only static states and indicates the amount of time spent in each of these static states as part of the 
transition. Specifically for this example, 

t 1 "? 1 = max(T ); t l ^ } = min(T .) - max(T ) ; and t 1 ~J J = t - min(T ). 
i v i->k k k -> j i -> k J cycle k -» j 

5 Hazards are introduced when they act to increase (in the case that we are seeking the maximum body 
voltage) or decrease (in the case that we are seeking the minimum body voltage) the steady-state 
body voltage that would result from the particular waveform being repeated indefinitely. 

This transition analysis determines the t J ~* k for the static states/. f 6 7_> * depends, 

5 in principle, on the number of switching events in a cycle which occur as a result of current flow 
1 Op through the target device. Each such switching event contributes an amount t swltch to t, J ^ k where 

tswitch is the approximate switching time of the FET. For the purposes of probabilistic body-voltage 
C estimation, it is assumed that if state 1 60 is accessible, exactly one switching event occurs per cycle. 
G In practice, the detailed results are not very sensitive to the number of assumed switching events or 
Li! the exact value of t switch .. For cases in which the on-state impact ionization current is large, any 

15* ^ k produces pinning at V r B ef = Y 6 zero yf° rward . For cases in which the on-state impact 
I ' 6 >0 

ionization does not dominate the current to the body in state 160, tf is usually much smaller than 
the tf associated with the static states. As a result, probabilistic body voltage has very little 
sensitivity to the exact value of tf . 
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From the t/^ k determined above, an effective amount of time (tf ) on the average 

per cycle that the FET is in the state / is determined by two sets of tf values, one to maximize and 
one to minimize the body voltage in accordance with equations 23 and 24: 

/ t eff ^max _ y pmax/ t k \max 
5 (t i } -^k^k^i > (23) 

{tf ) min = I j?k P^ in k (tf R ) min (24) 

Referring next to Fig. 7, the steps executed by a computer (not shown) executing 
1 (O software to implement the techniques described above will be described. In step 7 1 0, one or more 
m device models (n or p channel Field Effect Transistor models) are obtained. In step 720, the device 
£7 models are abstracted to generate one or more simplified electrical descriptions, each corresponding 
W to one of the transistors. As discussed above, the abstracting step may include abstracting each 
O device model to obtain displacement voltages d ; 72 1 , steady-state reference voltages Vi Zer0 722, and 
1 |*t forward reference voltages V t forward 723 , for each such model. Time constant characterizations 724 
2 are conditionally abstracted in the event that probabilistic determination is employed. 
*y In step 730, the circuit topology is checked to generate one or more sets of accessible 

□ states, with each set corresponding to one of the transistors and being indicative of whether a 
connection between a source or a drain of a transistor and either a power supply or ground exists. 
20 The conditional step 735 of performing active net tagging on each of net in the circuit to determine 
whether any of nets will switch with regular frequency is performed where modified accessibility 
analysis is employed. An alternative conditional step 736 of calculating signal probabilities and 
timing windows from the time constant characterizations by boolean and timing analysis is 
performed where probabilistic determination is employed. 
25 In step 740, sets of reference state body voltage minima and reference state body 

voltage maxima are determined for each of the transistors based on corresponding simplified 
electrical descriptions and corresponding sets of accessible states. As discussed above, this step may 
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include determining the sets of reference state body voltage minima and maxima based on 
corresponding displacement voltages d { and on corresponding set of accessible states 741, on 
corresponding displacement voltages d i? corresponding sets of accessible states, and from the results 
of the active net tagging 742, or on corresponding displacement voltages d i5 corresponding sets of 
5 accessible states, and from the results of the calculated signal probabilities and timing windows 743, 
depending on the type of analysis being performed. 

Finally, in step 750, the target state body voltage minima and a target state body 
voltage maxima are ascertained, one for each transistor, based on the determined sets of reference 
state body voltage minima and reference state body voltage maxima. Full uncertainty estimation 
10 751, accessibility analysis 752, modified accessibility analysis 753 or probabilistic analysis 754 
O may be employed. 

m In order to take advantage of the ability to accurately determine target state body 

J? voltage minima and maxima,, techniques for analyzing an electrical property of a digital partially 
fU depleted silicon-on-insulator circuit are now described. In one illustrative embodiment, the 
1 |fi technique involves ascertaining a target state body voltage minimum and a target state body voltage 
L maximum for each transistor in the circuit, establishing an initial condition for the circuit by 
5; selecting either the target state body voltage minimum or the target state body voltage minimum for 
5 each transistor, applying a voltage to the circuit, and measuring the sought after electrical property 
JSj of the circuit. Techniques for analyzing two exemplary electrical properties, switching delay and 
20 noise, are now discussed. 

First, the body- voltage estimates described above are applied to a static transistor- 
level timing analysis engine, or a SOI Static Timing Analyzer ("SOI-STA"), in order to analyze the 
electrical property of switching delay. The results for full-uncertainty,accessibility and probabilistic 
analysis are considered, although modified accessibility could also be easily applied in the context 
25 of static timing analysis. The design is partitioned into channel-connected components (CCCs) for 
analysis, as is traditionally done in static transistor-level tools. SOI-STA utilizes a breadth-first 
search (BFS) of the resulting timing graph, which ensures that the arrival times on all of the inputs 
are known when the delay through the CCC must be calculated. This enables detailed body-voltage 
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estimates on all of the FETs of the CCC to be used in delay simulation. Delay propagation through 
a CCC occurs as a result of a single switching event on an input (i. e., simultaneously switching 
inputs are not considered). Depending on the degree of knowledge of past switching activity (no 
knowledge, quiescent, steady switching, steady switching with knownsignalprobabilitiesand arrival 
5 times), one of the techniques discussed above (full-uncertainty, initial-condition, accessibility, 
detailed) is employed to determine the minimum and maximum possible value the reference body 
voltage can have for each FET of the CCC under analysis. For the given sensitization of the CCC 
for delay calculation, each FET is in a known state, the minimum or maximum body voltage of 
which can be determined by a displacement from the reference voltage. The body voltage values 
10 are then used as the "initial conditions" for the required delay simulation, 
p Early-mode calculation for rising transitions sensitizes the nFETs of the pull-up path 

m to be maximum and (to reduce the "fight" during switching) the nFETs of the pull-down path to be 
;D minimum; similar calculation for the pFETs of the pull-up path to be minimum and the pFETs of 
ftj, the pull-down path to be maximum. 
1 5g This same sensitization applies to late-mode fall transitions. Early-mode calculation 

!L. for falling transitions sensitizes the nFETs of the pull-down path to be maximum and the nFETs of 
ffl the pull-dup path to be minimum, and falling transitions sensitizes the pFETs of the pull-down path 
S to be minimum and the pFETs of the pull-dup path to be maximum. 

JSSRi. 

5l The same sensitization applies to late-mode rising transitions. The fanout CCCs are 

20 included in the delay calculation to improve the delay accuracy. To maximize device capacitance 
for late-mode, the maximum nFET and minimum pFET body voltages are used. To minimize device 
capacitance for early-mode, the minimum nFET and maximum pFET body voltages are used. 

SOI-STA propagates full piecewise-linear ("PWL") waveforms. One of the 
complexities of BFS timing analysis is determining the late- and early-mode waveforms at the output 
25 of each CCC. It can sometimes be the case that the waveform with the maximum (minimum) delay 
(as measured by the 50-percent point) is not the waveform with the slowest (fastest) slew (as 
measured by the 20-to-80-percent rise-fall time). In these cases, a "hybrid" waveform is propagated, 
chosen with the largest slew as the late-mode waveform and translate it in time so that is has the 
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largest delay, or with the smallest slew as the early-mode waveform and translate it in time so that 
is has the smallest delay. SOI-STA also propagates signal probabilities using assumptions of spatial 
and temporal independence, borrowing from similar techniques in static power analysis. If 
probabilistic body-voltage estimation is used, once the signal probabilities and arrival times are 
5 known at the inputs of a CCC ? these probabilities are translated into FET signal probabilities and 
arrival time values. 

Referring next to Figure 8, a 4-2 compressor circuit from a tree multiplier design is 
examined as an example of the advantageous use of the static transistor-level timing analysis 
presented herein, with the long delay path, I4-Ap-C-D-E-F-S, and the short delay path, I2-C-E-S, 
10 being considered. In the dynamic analysis, input waveforms of "50% duty cycle" are sued to 
rj sensitize these paths, corresponding to signal probabilities of 0.5 on all of the inputs in SOI-STA 
[£ detailed steady-state analysis. In Table III, the initial-condition delays determined by circuit 
C simulation are compared with the SOI-STA-determined full-uncertainty delays. 



TABLE III 




SOI-STA full-uncertainty 




SPICE initial 








Short Path 


Long Path 


Short Path 


Long Path 




rise 


fall 


rise 


fall 


rise 


fall 


rise 


fall 


12 


100 


100 


100 


100 


100 


100 


100 


100 


14 


100 


100 


100 


100 


100 


100 


100 


100 


Ap 






302 


491 






235 


315 


C 


301 


317 


866 


712 


326 


379 


476 


432 


D 






1057 


1079 






723 


553 


E 


526 


533 


1421 


1359 


1390 


762 


785 


943 


F 






1620 


1697 






1112 


936 


S 


672 


685 


1926 


1829 


1790 


837 


1104 


1174 
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Table IV does a similar comparison between the steady-state delays determined by 
circuit simulation (after more than 50,000 cycles of simulation) and the SOI-STA-determined 
probabilistic ("detailed")andaccessibilitysteady-state delays. In all cases, the SOI-STA-determined 
delays bound the SPICE delays. One should also notice the considerable reduction in uncertainty 
5 between the full-uncertainty and probabilistic steady-state delays, as the component of this 
uncertainty due to body voltage variation is noticeably reduced. 



10 


TABLE IV 




SOI-STA probabilistic 


SOI-STA accessibility 


SPICE steady-state 






Short Path 


Long Path 


Short Path 


Long Path 


Short Path 


Long ] 


Path 




12 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 . 


100 


100 




14 


100 


100 


100 


100 


100 


100 


100 


100 


100 \ 


100 , 


100 . 


100 ' 


m 


Ap 






267 


409 






296 


449 






266 


362 




C 


464 


485 


706 


605 


389 


365 


804 


673 


488 ; 


_503 ■ 


594 


576 




D 






897 


877 






1007 


999 






784 


794 




E 


970 


985 


1168 


1146 


733 


746 


1332 


1287 


1040 


1017 


1059 


1099 




F 






1363 


1409 






1540 


1590 






1246 


1352 


2ffl 


S 


1372 


1347 


1610 


1534 


1026 


986 


1812 


1737 


1384 


1365 


1508 


1483 



The body-voltage estimates described above are next applied to transistor-level static 
noise analysis, in order to analyze the noise immunity of the circuit. In this case modified 
accessibility analysis is considered for estimatics, although full uncertainty, accessibility, or 
25 probabilistic estimation could be applied to the circuit. As in the static time case, the approach 
involves decomposing the design into a collection of channel-connected components ("CCCs"), 
transistors that are connected together through their sources and drains. The maximum noise that 
is possible on each net is calculated as a time-domain waveshape. This worst-case noise considers 
all possible noise sources: leakage, charge-sharing noise, coupling though the interconnect, and 
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power-supply noise. This is done with a careful choice of vectors on the driving CCCs, referred to 
as the sensitization, which produces this worst case noise. 

Noise can also propagate from CCC-inputto CCC-output (propagated noise). Noise 
failures are determined by the noise stability, a type of AC noise margin analysis, of each CCC given 

5 the worst case noise appearing at its inputs. This involves calculating the transient sensitivity of the 
output noise with respect to the dc-level of the input noise. While the special considerations 
associated with applying static noise analysis to PD-SOI circuits is discussed herein, the article by 
K.L. Shepard et al., "Harmony: Static Noise Analysis for Deep-Submicron Digital Integrated 
Circuit" IEEE Trans. CAD 1 132-1 150 (1999), the disclosure of which is incorporated by reference 
1 0 herein, provides more details on static noise analysis generally. 

p There are two important considerations for static noise analysis for PD-SOI circuits. 

S The first is that all of the body voltages must be initialized as part of the CCC analysis used to 
calculate the noise on each node as well as the CCC analysis used to determine the noise stability 

flj of the gate. These body voltages are calculated using the modified accessibility analysis described 
1 !g above The approach here is straightforward and depends on the noise type being calculated. Two 

'1 types of noise on each CCC output, V H noise, which is noise that pulls the output down from the 

ffl supply level, and V L noise, which is noise that pulls the output up from ground, are calculated. 

5 When calculating V H noise on a CCC output (or when verifying noise stability in the 

;S! case that V H noise is introduced at the output), all of the devices in the pull-down stack are initialized 
20 to maximize the device strength by minimizing threshold voltages (maximum body voltages for 
nFETs and minimum body voltages for pFETs). This increases the strength of these devices in 
introducing noise. Similarly, all the devices in the pull-up paths are initialized to minimize device 
strengths by maximizing threshold voltages (maximum body voltages for pFETs and minimum body 
voltages for pFETs). This reduces the strength of these devices in maintaining the output at the logic 
25 high level. For V L noise, the situation is just the opposite. Devices in the pull-down path are 
weakened (minimum body voltages for nFETs and maximum body voltages for pFETs), while 
devices in the pull-up path are strengthened (maximum body voltages for nFETs and minimum body 
voltages for pFETs). 

NY02:2664023 Page -28- 



AP 32158 -070050.1280 

The second important consideration in the static noise analysis of PD-SOI circuits 
is that in calculating both propagated noise and in verifying noise stability, a parasitic bipolar 
leakage current may also be activated. Because of this, to achieve the worst-case, the arrival time 
of the input noise must be varied to maximize the output noise; that is, make sure that the noise 
5 introduced by parasitic bipolar leakage is superimposing maximally with the propagated noise 

Referring next to Figs. 9a and 9b, a simple circuit is examined as an example of the 
advantageous use of the noise analysis presented herein. Figure 9a illustrates a dynamic simulation 
of a circuit. Figure 9b shows the CCC decomposition used internally by the static noise analysis 
engine on the same circuit. In this partitioned circuit, the dynamic node is referred to as dyns and 
1 0 the output of the inverter as outs. 
O The static noise analysis results and dynamic simulation of the same circuit are 

m compared by vectors which correspond to the * ' worst-case" behavior detected in the static analysis. 
* In the circuit of Figure 4(a), the coupling noise introduced by the switching of an aggressor (G, H, 
rU or I) capacitively coupled to inputs A, C, and E, respectively, after a long period of switching history - 
1 fj is considered, called the prelude. Three different preludes are considered as shown in Table V below:/ - - 



TABLE V 


Prelude 






Node disposition 








A 


B 


C 


D 


E 


F 


1 


0 


0 


0 


0 


active 


0 


2 


0 


0 


active 


0 


0 


0 


3 


0 


0 


0 


0 


active 


active 



25 The clock period is lOnsec and the prelude is simulated for lOOmusec (so that the 

body has reached steady state) before the noise is applied. For the purposes of the dynamic 
simulation, a net that is noted as active in Table V is given a 0 - 1 transition at the middle of the 
evaluate cycle, every other evaluate cycle during the prelude. After the prelude, which is different 
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depending on the case considered, the "noise condition" shown in Table VI is applied, introducing 
V L coupling noise on node A for case 1, node B for case 2, and node C for cases 3 and 4: 



TABLE VI 


Case 


Prelude 


A 


B 


c 


D 


E 


F 


G 


H 


I 


1 


1 


0 


1 


0 


0 


0 


0 


fall 


1 


1 


2 


2 


0 


0 


0 


1 


0 


0 


1 


fall 


1 


3 


1 


0 


0 


0 


0 


0 


1 


1 


1 


fall 


4 


3 


0 


0 


0 


0 


0 


1 


1 


1 


fall 



This ^noise condition" is identical to the one established by static noise analysis to produce the 
worst-case V H noise on node dynd in each case, 
in The results of the dynamic circuit simulations are shown in Figures 1 0a through 1 3 a. 

l|i In each of these graphs, the dynamic node dynd, the dynamic gate output node outs, the V L noise 
Lj being injected on the gate, and the body voltage of the FET receiving the V L noise on its gate is 
a shown. The dynamic gate is given a * * full-latch" keeper so that the failure of the case is immediately 
O evi dent with the switching of this latch. Cases 1 and 3 show noise failures (as evidenced by the false 
□ switching of the latch), while cases 2 and 4 do not. In this dynamic case, the full prelude is 
2(tj simulated before the noise condition is applied. 
P Figure 10b through 13b show the noise on node dyns computed by static noise 

analysis. It should be noted that in the static analysis, the feedback path is broken and replaced by 
a tie. In this case, these is no "prelude simulation." Instead, the body voltages of each FET are 
initialized by the tool with the maximum or minimum value as determined by modified accessibility 
25 analysis. The inputs noted as active in Table V (in addition to the clock) are considered as active 
nets in this analysis. In Case 1, a significant parasitic bipolar current is also present; as a result, the 
tool adjusts the arrival time of the noise on A to produce the maximum peak noise at dyns. 

Figures 10c through 1 3c show the sensitivity ofthe domino-gate inverter output with 
respect to the noise calculated on dyns. Failures (as determined by sensitivity exceeding one in 
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magnitude) are correctly predicted for Cases 1 and 3. The noise is so large in Case 1 that the 
sensitivity has two " 'peaks" as the inverter is biased through the "high-gain" region and back again. 

The noise characteristics are most strongly affected by the body voltages on FETs 
Ml, M3, and M5. The diamonds on the left axis of Figure 10b, lib, 12b, and 13b indicate the 
5 "initial condition" body voltages used for FETs Ml, M3, M5, and M5, respectively. The actual 
(from dynamic simulation) and estimated maximum values (from the static noise analysis tool) for 
these body voltages are compared in Table VII. 



TABLE VII 


Transistor 


Actual Body 
voltage (state 5) 


Estimated max body 
voltage (state 5) 


M3 for case 2 


0.05 V 


0.255 V 


M5 for case 3 


0.720 V 


0.787 V 


M5 for case 4 


0.11 V 


0.255 V ' 


Ml for case 1 


0.95 V 


0.97 V 



J* In Case 1 , transistor Ml 3 charges the internal node between transistors M 1 and M2 

2 (JJ to V DD every precharge phase. This configuration would be common in bulk CMOS to mitigate 
y charge-sharing noise but is potentially disastrous for SOI since it can allow the body of M 1 to easily 
charge to a very high bias, nearly 1 V. When brought into state 150 by the noise condition, Ml 
draws bipolar leakage current. The high body voltage also reduces the threshold voltage, making 
the FET more sensitive to the coupling noise on the gate. 
25 Case 3 is similar, but the fact that net E is under steady switching means that both 

states 120 and 140 are accessible for FET M5. State 2 accessibility limits the maximum body 
voltage possible to only about 0.75 V. In this case, the body voltage is not high enough to result in 
significant parasitic bipolar leakage but the threshold voltage is still reduced enough to result in 
functional failure. 



NY02:266402.3 



Page -31- 



AP 32158 -070050.1280 

In Case 2, transistor Ml 4 discharges the internal node between transistors M3 and 
M4 every cycle. From the point of view of the modified accessibility analysis, this is equivalent to 
Case 4, in which both E and F are active. The conservative upper bound overestimates the body 
voltage. The main reason for this discrepancy is the neglect of state 160 which is exerting a 
downward pressure on the body voltage which is not (in the interests of conservatism) considered 
in the modified accessibility analysis. 

The foregoing merely illustrates the principles of the invention. Various 
modifications and alterations to the described embodiments will be apparent to those skilled in the 
art in view of the teachings herein. Gate leakage is an emerging influence on the body voltage and 
could easily be incorporated in to the estimation technique. It will thus be appreciated that those 
skilled in the art will be able to devise numerous systems and methods which, although not explicitly 
shown or described herein, embody the principles of the invention and are thus within the spirit and 
scope of the invention. 
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CLAIMS 

I claim: 



11. A method for statically estimating a body voltage of one or more transistors which form 

2 digital partially depleted silicon-on-insulatorcircuit having a predetermined circuit topology 

3 comprising said one or more transistors and one or more nets connecting said transistors, 

4 comprising the steps of: 

5 a. obtaining one or more device models, each corresponding to one of said one or more 

6 transistors; 

|g b. abstracting each of said device models to generate one or more simplified electrical 

g descriptions, each corresponding to one of said one or more transistors; 

$} c. checking said predetermined circuit topology to generate one or more sets of 

Iflj accessible states, each set corresponding to one of said one or more transistors -^and 

1 £; being indicative of whether under any allowable switching activity, the source, gate 

1 1^ or drain could be high or low. 

If! d. determining one or more sets of reference state body voltage minima and reference 

1 J state body voltage maxima, one for each of said one or more transistors, based on 

1 g corresponding simplified electrical descriptions and corresponding sets of accessible 

16 states; and 

17 e. ascertaining one or more target state body voltage minima and target state body 

18 voltage maxima, one for each of said one or more transistors, based on said 

19 determined sets of reference state body voltage minima and reference state body 

20 voltage maxima. 

1 2. The method of claim 1 , wherein said device models are selected from the group consisting 

2 of an n channel Field Effect Transistor model and a p channel Field Effect Transistor model. 
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1 3. The method of claim 2, wherein said abstracting step comprises abstracting each of said n 

2 channel Field Effect Transistor models, if any, and each of said p channel Field Effect 

3 Transistor models, if any, to obtain one or more displacement voltages d 1? and one or more 

4 steady-state reference voltages V, Zero , for each of said device models. 

1 4. The method of claim 3, wherein said corresponding simplified electrical descriptions 

2 comprise said displacement voltages d i; , and said determining step comprises determining 

3 said sets of reference state body voltage minima, and reference state body voltage maxima 

4 based on corresponding displacement voltages d, and corresponding set of accessible states. 

0 5. The method of claim 4, wherein said ascertaining step comprises ascertaining one or more 
M target state body voltage minima and one or more target state body voltage maxima, based 

Sir I 

; |P on said determined reference state body voltage minima and maxima using full uncertainty 

M estimation. 

;L 6. The method of claim 2, wherein said abstracting step comprises abstracting each of said n 
channel Field Effect Transistor models, if any, and each of said p channel Field Effect 

J Transistor models, if any, to obtain one or more displacement voltages d^, steady-state 

k: reference voltages Vj Zero , and forward bias reference voltages Vf°™ rd , for each of said device 

5 \\ 

5 models. 

1 7. The method of claim 6, wherein said corresponding simplified electrical descriptions 

2 comprise said displacement voltages d { , and said determining step comprises determining 

3 said sets of reference state body voltage minima, and reference state body voltage maxima 

4 based on corresponding displacement voltages d h and corresponding set of accessible states. 



1 

2 



8. The method of claim 7, wherein said ascertaining step comprises ascertaining one or more 
target state body voltage minima and one or more target state body voltage maxima, based 
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3 on said determined reference state body voltage minima and maxima using accessibility 

4 analysis. 

19. The method of claim 6, further comprising the step of performing active net tagging, after 

2 said checking step, on each of said one or more nets to determine whether any of said one 

3 or more nets will switch with regular frequency. 

1 10. The method of claim 9, wherein said corresponding simplified electrical descriptions 

2 comprise said displacement voltages d 1; , and said determining step comprises determining 

3 said sets of reference state body voltage minima, and reference state body voltage maxima 
|i based on corresponding displacement voltages d i<5 and corresponding set of accessible states, 
J . and from said determination from said active net tagging. 

|J 11. The method of claim 1 0, wherein said ascertaining step comprises ascertaining one or more 

S3 target state body voltage minima and one or more target state body voltage maxima, based 

% on said determined reference state body voltage minima and maxima using modified 

|t accessibility analysis. 

J 12. The method of claim 6, wherein said abstracting step further comprises abstracting each of 

2 said n channel Field Effect Transistor models, if any, and each of said p channel Field Effect 

3 Transistor models, if any, to obtain one or more corresponding time constant 

4 characterizations, one for each of said device models. 

1 13. The method of claim 1 2, further comprising the step of calculating signal probabilities and 

2 timing windows from said time constant characterizations, after said checking step, on each 

3 of said one or more nets, wherein said signal probabilities are determined by boolean 

4 analysis, and said timing windows are determined by timing analysis. 
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1 14. The method of claim 13, wherein said corresponding simplified electrical descriptions 

2 comprise said displacement voltages d i; , and said determining step comprises determining 

3 said sets of reference state body voltage minima, and reference state body voltage maxima 

4 based on correspondingdisplacement voltages d ij? corresponding sets of accessible states, and 

5 from said calculated signal probabilities and timing windows. 

1 15. The method of claim 14, wherein said ascertaining step comprises ascertaining one or more 

2 target state body voltage minima and one or more target state body voltage maxima, based 

3 on said determined reference state body voltage minima and maxima using probabilistic 

4 analysis. 

f 16. A method for analyzing an electrical property of a digital partially depleted silicon-on- 

jt 1 insulator circuit having a predetermined circuit topology comprising one or more transistors 

QJ and one or more nets, comprising the steps of: 

j£ a. ascertaining a target state body voltage minimum and a target state body voltage 

2L minimum for each of said transistors in said circuit; 

f b. establishing an initial condition for said circuit by selecting either said target state 

J body voltage minimum or said target state body voltage minimum for each of said 

% transistors in said circuit; 

9 c. applying a voltage to said circuit; and 

1 0 d. measuring said electrical property of said circuit. 

1 17. The method of claim 1 6, wherein said electrical property comprises a switching delay and 

2 said measuring step comprises measuring a delay between a switching input and a switching 

3 output as a constituent simulation for static timing analysis. 
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1 18. The method of claim 1 6, wherein said electrical property comprises noise and said measuring 

2 step comprises measuring noise on one or more nets in said circuit as a constituent 

3 simulation for static noise analysis. 

1 19. The method of claim 16, wherein said ascertaining step comprises ascertaining one or more 

2 target state body voltage minima and one or more target state body voltage maxima, using 

3 full uncertainty estimation. 

1 20 . The method of claim 1 6, wherein said ascertaining step comprises ascertaining one or more 

2 target state body voltage minima and one or more target state body voltage maxima, using 
gj accessibility analysis. 

f ■ 21. The method of claim 1 6, wherein said ascertaining step comprises ascertaining one or more 

!2J target state body voltage minima and one or more target state body voltage maxima, using 

p| modified accessibilityanalysis. 

CP 22. The method of claim 1 6, wherein said ascertaining step comprises ascertaining one or more 

3 target state body voltage minima, and one or more target state body voltage maxima, using 
3 probabilistic analysis. 
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ABSTRACT 



Techniques for estimating a body voltage of one or more transistors which form 
digital partially depleted silicon-on-insulator circuit, and for using estimated voltages to analyze 
electrical properties of the circuit, are disclosed. In one technique, device models are obtained and 
5 abstracted to generate simplified electrical descriptions of the transistors. The circuit topology is 
checked to generate sets of accessible states for the transistors that are indicative of whether a 
connection between a source or a drain of a transistor and either a power supply or ground exists. 
Next, sets of reference state body voltage minima and reference state body voltage maxima are 
determined for each of the transistors based on corresponding simplified electrical descriptions and 
Ol 0 corresponding sets of accessible states. Finally, a target state body voltage minima and a target state 
m body voltage maxima are ascertained, one for each transistor, based on the determined sets of 
f: reference state body voltage minima and reference state body voltage maxima. 
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